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Abstract 

Phase separated borosilicate glass samples were prepared by applying various heat 
treatments. Using selective chemical etching we performed AFM measurement on 
the phase separated glass surfaces. A quantitative roughness analysis allowed us to 
measure precisely the dependence of the characteristic size of the phase domains 
on heating time and temperature. The experimental measurements are very well 
described by the theoretically expected scaling laws. Interdiffusion coefficients and 
activation energy are estimated from this analysis and are consistent with literature 
data. 
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1 Introduction 



Within the last 20 years Atomic Force Microscopy (AFM) has become a stan- 
dard tool in glass science [1]. Beyond the direct roughness measurements, AFM 
has been applied to study the glass structure, corrosion, fracture, polishing and 
progressively more to characterize various (organic or inorganic) coatings de- 
posited on glasses. Increasingly precise measurements have been obtained in 
the last years allowing to reach atomic resolution [2,3]. Apart from detailed 
studies of the scale invariant character of the roughness of fracture surfaces [4] 
or the recent description of fused glass surfaces in terms of frozen capillary 
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waves [5], AFM roughness analysis is often restricted to the simple determina- 
tion of the standard deviation of the height fluctuations, i. e. the RMS rough- 
ness. This RMS roughness indicator which only gives information on the out 
of plane surface fluctuations is however a rather poor information as compared 
to the available information in the entire image. The characterization of in- 
plane features, e.g. the typical size of patterns on the surface requires the use 
of the height auto-correlation function. 



In this paper we develop a quantitative analysis of AFM roughness data 
to study the kinetics of phase separation coarsening stage in glasses. AFM 
measurement has been used to study phase separation but mostly to image 
the surface of thin polymer films [6,7,8,9,10,11]. The development of phase 
separation in glasses has been intensively studied experimentally in the last 
decades[12]. Various techniques have been used to follow this phenomenon : 
visible or X scattering[13,14,15,16,17], electronic microscopy[18,19], infrared 
or Raman spectroscopy[20, 21,22,23]... Most of the times, these techniques are 
used to test the occurrence of phase separation and to describe its morphol- 
ogy (droplet or interconnected). Few studies (generally based on electronic 
microscopy or X-ray scattering) give a quantitative estimate of the volume 
fraction of the two separated phases and sometimes of the typical size of the 
domains. X-ray and visible light scattering techniques have been intensively 
used to characterize the kinetics of phase separation, from the early stage of 
binodal or spinodal phase decomposition up to the scaling behavior of the 
coarsening stage[12]. More recently, coherent X-ray techniques have been used 
to probe the dynamics of fluctuations during domain coarsening of a sodium 
borosilicate glass [24]. 



In the present work we focus our study on a sodium-potassium borosilicate 
glass. Alkali borosilicate glasses [12] are known for their large composition do- 
main allowing phase separation. After heat treatment and etching we perform 
AFM surface measurements in order to reveal the internal structure. The paper 
is organized as follows. We first describe the experimental methods: melting, 
sample preparation, heat treatments, chemical etching and AFM measure- 
ments. We then present the topographic results obtained by AFM, giving a 
special attention to the statistical analysis. We finally compare our results with 
theoretical descriptions for phase separation kinetics and we extract estimates 
of interdiffusion coefficient and activation energy. 



2 



2 Experimental 



2. 1 Material 



A sodium potassium borosilicate glass of weight composition SiC>2 70%, B2O3 
25%, Na 2 2.5% and K 2 2.5% has been prepared as follows: i) The raw 
materials are mixed and melted at 1550°C in a 800 cm 3 platinum crucible; ii) 
The glass is then refined during 2 hours in order to obtain an homogeneous 
glass; iii) After quenching, glass is rapidly transferred to an electric furnace to 
be annealed during lh at 630° C, in order to relax the internal thermal stresses; 
iv) The glass is finally cooled to room temperature by opening the furnace 
door. After preparation, a cylindrical sample is core drilled and then cut into 
centimeter thick slices. At that stage, the glass is no longer transparent. Its 
opalescence is a signature of phase separation taking place over length scales 
comparable to the wavelengths of visible light. The major phase is expected 
to have a concentration close to pure silica while the minor phase concentrates 
the other constituents[12]. 



2.2 Thermal treatments 



Two types of thermal treatments have been performed. Four thick slices are 
first submitted to a thermal treatment at 650°C of duration lh, 4h, 16h and 
64h respectively. Square samples of 0.5mm thickness are cut from the slices 
and then mechanically polished to obtain mirror smooth surfaces. Following 
this first method, the result of bulk phase separation can be observed at the 
surface. A second sample preparation has consisted of inverting the sequence: 
first cutting the square samples from a glass slice, polishing them and only 
then performing the thermal treatments. Following this second method, we can 
observe the influence of a surface on phase separation. Two series of thermal 
treatments have been performed, one at 650°C with seven different durations 
(lh, 2h, 4h, 8h, 16h, 32h, 64h and 96h) and a second one at six different tem- 
peratures (625°C, 650°C, 660°C, 670°C, 675°C and 680°C) during 16h. The 
first series allows us to study the influence of thermal treatment duration on 
phase separation morphology. The second one allows investigation of the in- 
fluence of temperature. In all cases, the samples are heat-treated in an electric 
furnace and cooled to room temperature outside the furnace. 
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2.3 Etching 



A selective attack with a slightly acid treatment (pH 5) is performed. The 
borate rich phase is thus eliminated in surface due to its very low chemical 
durability while the silica rich phase remains almost unaffected. The resid- 
ual surface roughness then makes it possible to reveal the phase separation 
morphology as the composition contrast induces a height contrast. The alkali 
borate rich phase is thus expected to appear as depressions compared to the 
almost unmodified mirror surface corresponding to the silica rich phase. 

2.4 AFM measurements 

The surface morphology of all samples have been characterized by AFM height 
measurement in tapping mode (TM-AFM), using a Nanoscope III A from 
Digital Instruments with Al coated tip (BudgetSensors - model BS-Tap 300 
Al) . The Al coating thickness is 30nm, the resonant frequency is 300 kHz and 
the stiffness constant is 40N/m. The images have been recorded at a scan 
frequency between 0.8 and 1Hz for a resolution of 512 x 512 pixels. For each 
samples (i.e. for each thermal treatment conditions), we perform at least 6 
AFM images with 3 different scan areas. Depending on the size of the phase 
separation domains, the scan areas used are 1 x 1/im 2 , 2 x 2/im 2 , 4 x 4/im 2 
and 8 x 8fj,m 2 . 



3 Results 

3.1 Morphology of phase separation domains 

For all samples (cf. Fig. 1 and Fig. 2), we observed an interconnected mor- 
phology. On all these images, the grey scale is encoded such that the silica 
rich phase appears in clear colors and the borate rich phase appears in dark 
colors. On Fig. 1 (respectively on Fig. 2), one can easily see the influence of 
time (resp. temperature) of thermal treatment on the morphology: at a given 
temperature (resp. time), the characteristic size of spinodal phase separation 
domains increases with annealing time (resp. temperature). 

As clearly shown on Fig. 3, we obtained patterns exhibiting a striking self- 
similar character. In particular with a well chosen set of thermal treatment 
duration and AFM image size, we obtain almost indistinguishable images. 
This self similar character simply reflects that space and time can be rescaled 
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altogether. We show below that AFM measurements allow us to go beyond this 
qualitative picture. In particular we will evidence that the kinetics of the phase 
domains obeys the expected scaling law £ oc t 1 / 3 between the characteristic 
size £ and the duration t of the thermal treatment (at a given temperature) 
in the coarsening stage of the phase separation process[12]. 



3.2 Statistical analysis of AFM images 



In the following we use AFM measurements on etched surfaces to estimate the 
characteristic size of the phase separated domains. AFM measurements give 
an immediate access to the height field. To get information on the growth of 
the phase domains we will exploit the height autocorrelation function. 



3.2.1 Determination of the radial-mean auto- correlation function 

In order to evaluate the characteristic length of phase separation domains for 
all our samples, we first calculate for all TM-AFM height images the normal- 
ized auto-correlation function C(R), given by the following equation: 

r<tB\ - {h(x)-h(x + R)) 3 

C[R) - (h(x).h(x)) s (1) 



where h(x) is the height in nanometers at a point of coordinate x(r, 6) and 
with C(0) = I. On Fig 4, we give an example of the three dimensional 
representation of the center part (0.5x0.5/zm 2 ) of the auto-correlation function 
C (R) calculated from the height matrix extracted from an AFM height images 
(4 x 4/im 2 ). On this particular 3D representation and also for all other images, 
we obtain results characteristic of an isotropic random distribution of domains: 
a central peak surrounded by a few circular rings, corresponding to oscillations 
of vanishing amplitude. The isotropic character of the fluctuations allows for 
the use of an average auto-correlation function dependent on the only distance: 

g(r) = (C(r,9)) e (2) 



Fig. 5 shows the auto-correlation function for a heat treatment of 16h at 650°C 
performed on two series of samples where the heat treament has been applied 
before or after cutting and polishing operations. The evolution of g(r) as a 
function of the duration of heat treatment at 650° C and respectively as a 
function of temperature for a 16h heat treatment is given on Fig. 6 and Fig. 7 
respectively. On these two figures, g(r) corresponds to the mean value of the 
auto-correlation function of the six AFM images associated with one sample. 
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3.2.2 Determination of two correlation lengths 

The second step of our analysis of AFM height data consists of extracting 
for each sample a characteristic length from auto correlation function g(r). 
We discuss here two typical lengths, chosen for their physical meaning with 
respect to phase separation phenomena. The first one, called £o.5> is the value 
of r at g(r) = 0.5 (#(£0.5) — 0.5)) and corresponds to the mean size of the 
phase separation domains. The second one, called £, m in, is the value of r when 
g(r) reached its first minimum (g(£ m i n ) = min(g(r))) and corresponds to the 
mean distance between two domains of different composition. The evolution of 
these characteristic lengths with heat treatment duration and temperature are 
summarized in tables 1 and 2 and displayed in Fig. 9 and Fig. 10 respectively. 
These results are discussed in next section in the light of theoretical predictions 
for the kinetics of the coarsening stage of phase separation in glasses. 



4 Discussion 



The process of phase separation in a spinodal region is usually divided into 
two stages. The first one corresponds to the phase decomposition into domains 
at the two limit concentrations. The second one is the coarsening process, 
corresponding to the progressive growth of the domains, driven by surface 
tension and limited by diffusion. During this coarsening regime, the structure 
remains statistically self-similar and the evolution of the morphology is entirely 
controlled by the growth of the characteristic length of the domains which is 
expected to follow the Lifshitz-Slyozov- Wagner scaling: 

R^{Ktf/\ K*^, (3) 



where, within the classical mean field theory of coarsening[12,25] 7 is the 
interfacial tension, D the interdiffusion coefficient and v a molecular volume. 
This result applies in case of conservation of the order parameter. Following 
Huse[26], it can be recovered via a simple scaling argument. Let us consider 
an interface of radius of curvature R separating two domains with a volume 
concentration difference Ac. Under local equilibrium, the Laplace pressure 
■y/R balances the osmotic pressure /iAc (where \i is the chemical potential) 
so that /i ~ '-f/RAc. Assuming that fluctuations of the chemical potential 
are of the same order of magnitude as the chemical potential itself we have 
V/i ~ r y/R 2 Ac. The flux of particles can thus be written j = McVfi where 
M = D/kT is the particle mobility. Mass conservation finally gives for the 
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domain growth law: 



dR 
~dt 



J 



c 



R 2 AckT 



7D 



(4) 



from which we recover directly the scaling law (3). 

In the following we shall compare our experimental results within this frame- 
work. The different durations of heat treatments allow us to study the tem- 
poral scaling of the coarsening regime and beyond, to estimate the diffusion 
constant. In a similar way, data from different temperatures of heat treatment, 
combined with the Arrhenian scaling of the mobility, M = MQexp(—E in /kT) 
provides a rough estimate of the corresponding activation energy. 

4-1 Effect of surface 

We summarize on Fig. 5 the height autocorrelation functions obtained after 
a heat treatment of 16 hours at 650°C Two types of preparation were used. 
As described above, the heat treatment was applied either before (series A) 
or after (series B) cutting and polishing the samples. This allows to test the 
effect of a surface on the coarsening process. As shown on Fig. 5 the correlation 
functions obtained on the two series of samples are almost identical. The 
presence of a surface thus does not seem to influence the coarsening and surface 
measurements can be considered as representative of the bulk morphology. The 
results presented in the following have been obtained by statistical averaging 
over the entire set of samples. 

4-2 Effect of annealing time 

Let us note first that all samples described here were submitted to a pri- 
mary annealing thermal treatment of lh at 630° C. The latter is intended to 
relax thermal stresses but also induces phase separation in the present case. 
Additional thermal treatments thus only tend to coarsen the existing phase 
separated domains. As an illustration, we plotted on Fig 8 the auto-correlation 
functions g(r) of all samples after rescaling r/£ both for £ = £ .5 an d £ = £mm- 
As expected in the coarsening stage, we observe that all correlation functions 
fall onto a master curve, indicating that, regarding the morphology of phase 
separation, the only changing parameter is a characteristic length scale. 

We report on Fig. 9 the time evolution of both correlation lengths £0.5 and £ m j„ 
in function of the duration of thermal treatment t at temperature 650°C. In 
this plot, we use the cubic root of time for the abscissa in order to reveal the 
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expected scaling. In these coordinates we get a very good linear behavior over 
the entire time range experimentally studied (from lh to 96h). The time evolu- 
tion of the size of phase separated domains evaluated by AFM thus obeys very 
well the predicted scaling for coarsening. Extracting the leading coefficients 
by performing linear fits we obtain 7D ~ 1CT 18 where we took v ~ 10~ 29 m 3 
for a molecular volume. Assuming for the surface tension a value of order 
7 ~ lCT 2 J.m~ 2 this gives us for the value of the interdiffusion coefficient at 
650° D ~ 10 _16 cm.s _1 , which is consistent with literature [12]. 

4-3 Effect of annealing temperature 

The same analysis is repeated to analyze the evolution of the correlation 
lengths £0.5 and £ mi „ for 16h long heat treatments at temperatures varying 
from 625°C to 680°C We now test the Arrhenius dependence of the interdif- 
fusion coefficient. The results are reported on Fig 10 in coordinates log£ vs 
1/T. We obtain again a reasonable linear behavior indicating that the Arrhe- 
nius scaling gives a good approximation of the dependence of the interdiffusion 
coefficient with temperature. 

The values of c and d parameters obtained via the linear fit In £ = c — d/T 
are given in table 3 with their associated standard deviations both for £ = £ .5 
and for £ = £ m j n . Using d = E in /3RT (the factor 3 coming from the expo- 
nent 1/3 in Eq. (3)) we can give a rough estimate of the activation energy 
Ei n = 135 ± 45kJ.mol~ 1 . This value seems to be consistent with typical re- 
sults from literature: by electrical resistivity measurements in ion exchanged 
Corning 7740 borosilicate glasses, Garfinkel[27] obtained activation energies 
in the range 80 — 120kJ.mol _1 ; activation energies measured for the diffusion 
of sodium in silicate networks are also in the same range [28]. 



5 Conclusion 

Performing selective etching on the surface of phase separated alkali borosili- 
cate glasses, we revealed the phase separated domain by AFM. A quantitative 
statistical analysis of AFM data allowed us to extract the characteristic corre- 
lation lengths of the coarsening stage of phase separation. A systematic study 
of the dependence on time and temperature of the phase separated domains 
was performed. We recover the expected cubic root time scaling for the domain 
size temporal evolution and we could extract an estimate of the interdiffusivity 
coefficient. The dependence of the latter on temperature is consistent with an 
Arrhenius behavior, allowing to give a rough estimate of the activation energy 
for interdiffusion. The numerical vaues obtained for these two parameters are 
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moreover consistent with literature data. This study performed on a very sim- 
ple phase separated glass thus shows that, beyond its performance in imaging, 
AFM analysis can be used in a more quantitative way to characterize phase 
separation, and be considered as an alternative or complementary technique 
to visible or X-ray scattering traditionally used in that purpose. 
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Duration 


fo.5 (nm) 


£, min (nm) 


1 hour 


27 


97 


4 hours 


31 


121 


8 hours 


43 


144 


16 hours 


58 


195 


32 hours 


66 


238 


64 hours 


94 


332 


96 hours 


105 


422 



Table 1 

Correlation length (£0.5 and £ m i n ) vs duration of thermal treatment at 650°C. 



Temperature 


£0.5 (nm) 


t, min (nm) 


625°C 


47 


148 


650°C 


58 


195 


660°C 


51 


160 


670°C 


62 


191 


675°C 


63 


207 


680°C 


66 


230 



Table 2 

Correlation length (£0.5 and £ m j n ) vs temperature for 16 hours thermal treatment. 



ln£ = c- # 


C 


d 


£0.5 

£.min 


-11.34 ± 1.54 
-9.43 ±2.13 


(4.9 ±1.4) • 10 3 
(5.6 ± 1.9) • 10 3 



Table 3 

Parameters from linear curve fitting of In £0.5 and ln£ m j n as a function of the inverse 
of absolute temperatures (1/T) for 16 hours thermal treatments. 
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Fig. 1. AFM height images of borosilicate glass samples after thermal treatment at 
650° C during six different durations. 
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|e) 675 C-16h 



Fig. 2. AFM height images of borosilic. 
treatments at six different temperatures. 
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glass samples after 16 hours thermal 
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Fig. 3. AFM height images of three borosilicate glass samples. These images il- 
lustrate the similarity of morphology of the phase separation domains for various 
observation scales according to the duration of thermal treatment at 650° C. 
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Fig. 4. AFM height images (4 x of a borosilicate glass sample after 64 hours 

of thermal treatment at 650° C (left) and the associated three dimensional repre- 
sentation of the center part (0.5 x 0.5fim 2 ) of the auto-correlation function C(x) 
(right). 
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Fig. 5. Radial mean value of the auto-correlation function g(r) vs distance r (in nm) 
for borosilicate glass samples annealed 16 hours at 650°C. Samples A have been cut 
and polished after annealing treatment and conversely samples B have been cut and 
polished before annealing treatment 
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Fig. 6. Radial mean value of the auto-correlation function g(r) vs distance r (in nm) 
for borosilicate glass samples annealed during different durations at 650°C. 
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Fig. 7. Radial mean value of the auto-correlation function g(r) as a function of 
distance r (in nm) for borosilicate glass samples annealed during 16h at different 
temperatures. 
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Fig. 8. Radial mean value of the auto-correlation function (g(r)) as a function of 
the normalized distance for borosilicate glass samples after thermal treat- 

ment during different durations at 650° C : (left) r is normalized with the cor- 
relation length £o.5 (<?(£o.5) = 0.5) and (right) with the correlation length £ m j n 
{g{£,min) = min(g(r))). 
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Fig. 9. Evolution of the two correlation length (£0.5 and £ m j n ) w.r.t. the duration of 
thermal treatment at 650° C. 
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Fig. 10. Evolution of the two correlation lengths (£0.5 and £ r , 
ature for 16 hours thermal treatments. 
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